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I.  INTRODUCTION 


The  Naval  Research  Laboratory  (NRL)  Limited  Area 
Dynamical  Weather  Prediction  Model  is  an  evolving  research 
model  for  testing  new  methods  of  modeling  various  mesoscale 
phenomena.  The  model  has  shown  success  in  resolving 
topographic  and  coastal  features.  Also,  high  vertical 
resolution  in  the  lower  levels  of  the  model  has  allowed 
close  examination  of  processes  occurring  within  the  boundary 
layer.  Previous  studies  (Grandau  1992  and  Stewart  1992) 
have  demonstrated  the  NRL  model's  effectiveness  in  resolving 
mesoscale  and  boundary  layer  features  in  this  region. 

The  United  States  West  Coast  poses  a  big  modeling 
challenge  due  to  its  topographic  features  and  coastal 
mesoscale  phenomena  such  as  the  Catalina  eddy  and  the 
southerly  surge.  In  addition,  much  of  the  California  coast 
experiences  frequent  stratus  clouds  which  can  significantly 
affect  the  weather  in  the  region. 

The  stratiform  condensation  (or  non-convect ive)  process 
is  important  in  modeling  coastal  mesoscale  phenomena.  A 
very  simplistic  approach  is  to  represent  cloud  cover  as 
being  either  total  or  zero  for  each  grid  box  depending  on 
whether  gridpoint  relative  humidity  has  reached  100%  or  not. 
A  more  advanced  technique  would  be  to  parameterize  the 


process  on  the  subgrid  scale. 


Cloud  cover  (or  cloud 


fraction)  could  be  determined  more  realistically  as  a 
function  of  gndpoint  relative  humidity  and  a  relative 
humidity  threshold  or  critical  value. 

This  study  examines  the  incorporation  of  the  stratiform 
condensation  process  of  Sundqvist  et  al.  (1989)  into  the  NRL 
limited  area  mesoscale  model.  Previous  NRL  model 
simulations  by  Stewart  (1992)  included  cloud  fractions  as 
determined  by  the  method  of  Slingo  and  Ritter  (1985) .  These 
simulations,  hereafter  referred  to  as  the  "control"  case, 
will  be  compared  to  the  NRL  model  incorporating  the 
Sundqvist  et  al.  explicit  non-convect i ve  cloud  liquid  water 
scheme,  hereafter  referred  to  as  the  "test"  case.  The 
inclusion  of  cloud  water  as  a  variable  in  predicting 
stratiform  clouds  is  a  major  feature  of  Sundqvist 's  scheme. 

Section  II  describes  the  NRL  mesoscale  model  and  the 
extent  of  the  geographic  region  involved.  Section  III 
outlines  the  Sundqvist  et  al.  stratiform  parameterization 
technique.  Section  IV  is  the  regional  weather  scenario  for 
the  time  period  including  specific  mesoscale  phenomena. 
Section  V  describes  the  mesoscale  structure  of  the  control 
model  simulations.  Section  VI  evaluates  the  test  model  runs 
in  comparison  to  the  control  model,  and  Section  VII  includes 
conclusions  and  recommendations. 


2 


I  I  .  MODEL  DESCRIPTION 


The  Naval  Research  Laboratory  (NRL)  regional  weather 
prediction  model  is  a  quas i -hydrost at  1 c .  baroclinic, 
mesoscale  model  which  incorporates  cumulus.  boundary  layer, 
and  radiation  parameter izat ions .  This  limited  area  model  is 
most  appropriate  where  near-gradient  balance  of  large  scale 
motions  in  the  lower  troposphere  exists. 

Since  the  isobaric  coordinate  system  does  not  easily 
handle  varying  topography,  the  vertical  coordinate  sigma  (a) 
is  used.  This  is  defined  as  the  ratio  of  pressure  to  the 
surface  pressure. 

<X=(-£-)  (2.1) 

P# 

Specifics  of  the  NRL  model  are  detailed  by  Madala  et  al. 
(1987) .  The  basics  are  outlined  in  the  following  pages. 

A.  GRID- 

The  horizontal  grid  is  a  staggered  Arakawa  C~grid.  'This 
type  of  network  is  best  for  the  simulation  of  wind  field 
geostrophic  adjustment  and  conservation  of  integral 
properties.  General  curvilinear  horizontal  coordinates  with 
user  specified  horizontal  grid  spacing  is  used.  For  an  MX 
N  field  (i=l,2,...M;  j=l,2,...N),  temperature  (T) . 
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geopotential  height  (<$),  specific  humidity  (qJ.  and  sigma 


(a)  are  computed  at  mass  points  ii.j;,  with  u-velocity 
(east-west)  and  v-velocity  ( north -south)  computed  at  the 
midpoints  along  the  and  y-axis  respectively  )See  Figure 
1)  . 

The  horizontal  domain  is  a  103  X  91  grid  with  1/6  degree 
resolution  in  latitude  and  longitude  from  2E7  N  to  43°  N  and 
113°  W  to  13(7  W.  Figure  2A  shows  the  extent  of  this  region 
along  with  points  of  interest  referred  to  for  profile  and 
cross-sectional  plots  in  this  study  and  Figure  2B  displays 
main  geographic  points  of  reference.  Mode!  simulations 
employ  23  vertical  sigma  levels  as  shown  in  Table  1. 
Thirteen  of  these  layers  are  below  850  mb  ensuring  a  high 
vertical  resolution  within  the  boundary  layer. 
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Figure  1 

Horizontal  and  vertical  grid  network 
utilized  In  the  NRL  meeoBcale  model 

(Nodal*  et  al.  13S7) 
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US  120 


Figure  2B 

Significant  geographic  points. 
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TABLE  1 

MODEL  SIGMA  LEVELS 


Model  Level 


Slgaa  (0) 


B.  EQUATIONS- 


The  governing  primitive  equations  are  formulated  in 
surface  pressure  flux  forms  (i.e.  pB u .  psv.  etc.).  The 
dynamic  system  consists  of  seven  equations.  five  of  which 
are  prognostic,  and  two  which  are  diagnostic.  (See  Madala. 
198?  for  the  complete  form  of  the  equations) . 


u-  &  v-momentum  equations: 


{P»U) 


(2.2) 


(2.3) 


thermodynamic  equation: 


(2.4) 


moisture  continuity  equation: 


-ik  <JVtf 


(2.5: 


surface  pressure  tendency  equation: 


£<p.> 


(2.6) 
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hydrostatic  equation: 


continuity  equations: 


(2.71 


(2.8) 


A  closed  system  is  formed  for  the  dependent  variables  u,  v, 
T,  q,  p8  ,  ‘i .  and  d  (vertical  velocity)  . 


C.  TIME  I NTEGRATION- 


The  split-explicit  method  is  utilized  which  effectively 
splits  terms  in  the  prognostic  equations  into  two  parts: 
those  governing  the  Rossby  modes  and  those  governing  the 
faster  gravity  modes.  For  quasi-1 inear  gravity  modes.  the 
pressure  gradient  and  divergence  terms  vary  faster  than  the 
remaining  terms.  This  allows  that  part  of  the  equation  with 
these  remaining  terms  to  use  a  larger  time  step.  The  split 
equations  are  integrated  with  time  steps  for  their 


respective  CFL  criteria. 
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D.  PARAMETER I ZATIONS- 


The  model  parameterizes  a  number  of  physical  processes 
including  cumulus  precipitation,  planetary  boundary  layer 
(PBL)  processes,  and  radiation. 

Cumulus  parameterization  is  from  the  modified  Kuo  scheme 
(Kuo  1974) .  Surface  layer  parameterization  is  based  on 
Mom n-Obukhov  similarity  theory.  PBL  parameterization  is 
with  turbulent  kinetic  energy  (TKE)  closure  described  in 
Holt  and  Raman  (1988) . 

The  radiation  parameterization  incorporated  in  the  model 
is  the  Harshvardhan  et  al.  (1987)  scheme.  Stewart's  (1992) 
study  parameterizing  longwave  and  shortwave  effects  showed 
the  improvements  in  the  model's  ability  to  simulate  diurnal 
and  cloud-related  radiative  processes. 

The  cloud  parameterization  scheme  incorporated  in  the 
con  . o 1  model  by  Stewart  (1992)  is  based  on  cloud  fractions 
using  a  modified  Slingo  and  Ritter  (1985)  method  where 
average  layer  relative  humidities  are  compared  to  critical 
relative  humidity  values.  This  method  produces  stable  and 
convective  cloud  fractions  for  a  horizontal  grid  at  each 
sigma  layer.  The  clouds  are  diagnosed  as  either  stratiform 
or  cumulus. 

Due  to  the  frequent  presence  of  stratiform  clouds  along 
the  west  coast,  the  need  for  a  more  sophisticated  non- 
convective  parameterization  scheme  is  clear.  This 


stratiform  parameterization  is  discussed  in  detail  in 
Section  III. 

E.  INPUT  DATA- 

Data  for  the  period  0000  UTC  02  May  -  1200  UTC  03  May 
1990  was  taken  from  the  Navy  Operational  Global  Atmospheric 
Prediction  System  (NOGAPS).  Model  initialization  data  was 
retrieved  from  archived  Fleet  Numerical  Oceanography 
Center's  2.5  degree  global  analyses  and  horizontally  and 
vertically  interpolated  to  the  NRL  model  resolution  (Grandau 
1992) .  Fields  include  u-  and  v-  components  of  velocity, 
temperature  (T) ,  vapor  pressure,  sea  level  pressure,  and  sea 


surface  temperature  (SST) . 


III.  STRATIFORM  PARAMETERIZATION 


A.  DESCRIPTION- 

Increased  development  of  mesoscale  atmospheric  modeling 
has  resulted  in  a  more  involved  look  at  the  proper  treatment 
of  the  stratiform  condensation  process.  The  high  occurrence 
of  stratus  clouds  along  the  U.S.  West  Coast  and  its  effect 
on  regional  weather  patterns  makes  this  region  ideal  for 
this  type  of  study. 

One  assumption  typically  made  in  representing  stratiform 
condensation  is  that  the  gridpoint  relative  humidity  must 
reach  100%  in  order  for  condensation  to  occur.  This  results 
in  a  simplistic  approach  of  cloud  cover  either  being 
represented  as  zero  or  100%  at  individual  grid  points.  A 
more  realistic  treatment  of  this  case  would  be  a  subgrid 
scale  method  requiring  a  parameterization  of  the  process. 

Sundqvist  et  al.  (1989)  devised  a  treatment  of 
condensation  cloud  processes  for  convective  and  non- 
convective  precipitation.  This  study  will  only  incorporate 
Sundqvist 's  non-convect ive  scheme.  The  convective  case 
utilizes  the  scheme  by  Kuo  (1974)  adopted  to  account  for  the 
inclusion  of  cloud  water  as  a  prognostic  variable.  The 
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first  step  in  the  parameterization  of  convective  and  non- 


convective  precipitation  is  a  stability  check  of  the  grid 
column.  The  criterion  here  is  that  an  air  parcel  at  the 
surface  should  be  positively  buoyant  after  reaching  the 
lifting  condensation  level  (LCL) .  If  the  column  is 
conditional ly  unstable,  the  Kuo  scheme  (without  Sundqvist's 
modifications)  is  used.  If  not.  the  possibility  of 
stratiform  condensation  is  investigated.  The  Sundqvist  et 
al.  stratiform  condensation  treatment  is  described  in  the 
following  pages. 

B.  STRATIFORM  CONDENSATI 0N- 


For  stratiform  condensation  to  take  place,  a  relative 
humidity  threshold  value  (of  less  than  unity)  within  a  grid 
box  must  be  exceeded.  Parameterization  of  the  stratiform 
condensation  process  is  a  function  of  quantities  such  as 
stability,  cloudiness,  altitude,  and  type  of  surface.  The 
prognostic  equations  used  are  those  for  cloud  water  mixing 
ratio  (m) ,  temperature  (T) ,  and  specific  humidity  (q) . 

Changes  in  cloud  water  are  due  to  local  changes  in  m 
from  advection  and  diffusion  (A  .  )  ,  latent  heat  release  (Q)  , 
local  rate  of  release  of  precipitation  (P) ,  and  evaporation 
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of  cloud  water  to  vapor  advected  to  a  grid  box  where  no 


condensation  is  occurring  (Ec) .  In  equation  form. 


dm 

dt 


-Aa+Q-P-Ec 


(3.1) 


where  dm/dt  is  total  cloud  water  mixing  ratio  tendency. 

The  local  rate  of  release  of  precipitation  (P)  is  dependent 
upon  a  cloud  liquid  water  threshold  value  (m  r)  .  Efficient 
release  of  precipitation  occurs  when  m  exceeds  the  cloud 
fraction  multiplied  by  this  threshold  value. 


P-Cf/nll-e 


(3.2) 


c,  is  characteristic  conversion  rate  of  cloud 
particles  to  precipitation  size,  and 
b  is  the  cloud  fraction. 


The  cloud  fraction  is 
humidity  (RH) ,  a  relative 
condensation  (RHC),  and  the 
of  the  cloudy  part  (RHB) . 
value  of  1 )  . 


determined  from  the  relative 
humidity  threshold  value  for 
weighted  average  of  the  humidity 
(RH9  is  defined  to  equal  the 


i>=l- 


{RHa-RM 
N  (RH.-RHC) 


(3.3) 
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Typically,  RH c  values  are  empirically  derived.  Sundqvist 
et  al.  (1989)  used  RH  threshold  values  above  the  boundary 
1 ayer  of : 

0.75  —  over  land,  and 
0.85  —  over  the  ocean. 

In  the  boundary  layer.  RHc  is  assumed  to  linearly 
approach  unity  at  the  surface.  Subgrid-sca le  topography 
effects  are  taken  into  account  by  assuming  the  value  of  RHc 
to  be  0.1  lower  over  the  land  than  over  the  ocean.  Also,  to 
prevent  the  unrealistically  early  formation  of  cirrus 
clouds,  RHC  is  increased  asymptotically  toward  unity  for 
temperatures  less  than  238°  K.  Figure  3  depicts  the 
Sundqvist  et  al.  relative  humidity  threshold  profile. 

This  is  contrasted  with  the  approach  by  Slingo  and 
Ritter  (1985)  and  used  by  Stewart  (1992)  in  the  control 
model  simulations.  Here.  different  critical  relative 
humidity  values  are  computed  as  a  function  of  sigma.  For 
each  model  sigma  level,  a  critical  relative  humidity  is 
determined  from  the  equation: 

J2Hc=l+2  (<72-<j)  +y^T<7  (l-3o+2a2)  (3.4) 

In  order  to  account  for  the  lack  of  model -simulated 
clouds  near  the  coastal  regions,  Stewart  used  satellite 
imagery  to  propose  a  modification  of  critical  relative 
humidity  values  for  this  particular  2-3  May  1990  case  study. 
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Figure  4  shows  the  representations  of  critical  relative 
humidity  values  by  both  Slingo  et  al.  and  Stewart. 

For  stratiform  clouds,  the  evaporation  of  precipitating 
water  to  vapor  as  it  falls  through  subsaturated  air  iEr) 
depends  on  cloud  fraction  (b).  relative  humidity  (RH).  and  a 
layer-averaged  precipitation  rate  CP) . 

ET=kB(RHB-RH)  (l-b)jp  (3  5) 

kt  is  a  coefficient  expressing  the  instantaneous  evaporation 
of  advected  m. 

Temperature  changes  are  due  to  the  temperature  tendency 
from  advection.  diffusion.  and  radiation  C A T )  .  latent  heat 
release  (Q) ,  evaporation  of  precipitating  water  to  vapor 
(Er  )  .  and  evaporation  of  advected  cloud  liquid  water  (Ec). 

jZ=AT+(-±-)Q--£-(Ec+Ez)  (3.6) 

cp  in  this  equation  is  the  specific  heat  of  dry  air  at 
constant  pressure,  and  L  is  latent  heat  of  vaporization. 

The  change  in  specific  humidity  idq/dt)  comes  from  the 
effects  of  advection  and  diffusion  (A,),  latent  heat  release 
(Q) ,  and  evaporation  (Ec+Er). 

jZ=A9-Q+Ec+ET  (3.7) 
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Sigma 


Figure  4 

Critical  relative  humidity  profiles. 

(Slingo  &  Ritter  —  solid) 

(Stewart  —  dashed) 


(Stewart  1992) 


Combining  the  prognostic  equations  tor  temperature  and 
speed  ic  humidity  with  the  C 1  aus  1  us-C  1  apeyron  equation 
(Sundqvist  1988)  gives  an  expression  tor  Q, 


0= 


1  + 


(RH)  eL2 qa 
RCpT3 


(.3.8] 


where  M  is  the  convergence  of  available  latent  heat  given 


as  : 


wags.  a, 

g  RT2  P  dt 


(3.9) 


q.  is  saturation  specific  humidity, 

6  is  the  ratio  of  molecular  weight  of  vapor 
to  molecular  weight  of  dry  air  (=0.622),  and 
R  is  the  gas  constant  for  dry  air. 


Hence,  positive  advection  of  A,  and  pressure  tendency 
(dp/dt)  would  tend  to  increase  available  latent  heat 
convergence,  while  warm  advection  (A  T  >  0)  would  tend  to 

decrease  convergence . 

To  close  the  system,  the  tendency  of  RH  can  be  expressed 

by: 

d(RH)  _  2(1-1?)  ( RH„-RHC )  [(l-b)M+Er) 
dt  2g,(l-i?)  (RH,-RHc)  +  (2)  (3.10) 

At  this  point,  temperature  and  specific  humidity 

tendencies  (ST/dt  and  dq/9t)  can  be  computed.  The  mixing 
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ratio  tendency  ( dm/dt )  is  obtained  by  semi -imp  1 1 c 1 1  time 
integration  of  equation  (3.1),  which  may  be  rewritten  as: 


dm 

dt 


=RHS-P 


(3.11) 


-j^=RHS-cam[l-e  '  ***' 


(3.12) 


(RHS  =  right  hand  side) 


And  in  finite  difference  form: 


mD+1=mn-t+2At(RHS)  -2Atcja[l-e  ] 


(3.13) 


£ 


(3.14) 


Resulting  in  a  non-linear  equation. 


jc[l+Afc0(l-e-)ta)]--^'1--i:^S)At 


bmr 


(3.14) 


(3.15) 


which  may  be  solved  using  the  Newton-Raphson 
method . 

Table  2  provides  a  summary  of  various 
condensation  processes  described  by  the  tendency 
and  their  influences.  Water  vapor  condensing 


iteration 

stratiform 
equations 
into  cloud 
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TABLE  2 

SUMMARY  OF  STRATIFORM  CONDENSATION  PROCESSES 


So 


dt  * 


=  A  ♦  _Lq  - 

dt  T  °D 


\TP!titl1iW?5T3*TVT7tTtctrTT7T!WTI?fIfJ!7!t 


PROCESS 


SIONIPIOANT 

TERM 


q  condenses 
to  m 

Q 

m  advected  & 
evaporates  to  q 

precipitation 
evaporates 
to  q 

E, 

m  precipitates 
out 

P 

4.T.M  TENDENCY 


ql  T|  m| 


q|  T  |  m | 


q|  T  |  m- 


liquid  water  results  in  a  decrease  in  specific  humidity  (q). 


and  an  increase  in  temperature  (T)  and  cloud  liquid  water 
(ml .  When  cloud  water  is  advected  and  evaporated  into  water 
vapor,  specific  humidity  is  increased,  while  temperature  and 
cloud  liquid  water  are  decreased.  For  the  case  where 
precipitation  evaporates  to  vapor,  specific  humidity  is 
increased  and  temperature  is  decreased.  Finally,  for  cloud 
water  precipitating  out,  cloud  liquid  water  decreases,  with 
specific  humidity  and  temperature  remaining  constant. 
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IV.  WEATHER  SCENARIO 

A.  SYNOPTIC  SITUATION: 

1 .  Upper  Air- 

The  period  of  interest  for  this  study  is  0000  UTC  02 
May  1990  to  1200  UTC  03  May  1990.  A  high  pressure  ridge 
dominated  the  eastern  Pacific  region  and  slowly  intensified 
during  the  period.  A  closed  upper  level  low  that  was 
situated  over  southern  Arizona  initially  deepened  and 
subsequently  filled  and  moved  eastward  over  the  latter  part 
of  the  period. 

There  was  also  evidence  of  a  weak  shortwave  moving 
over  the  Pacific  Northwest  at  the  start  of  the  period. 
Figures  5A  to  5D  show  the  500  mb  pattern  for  the  period. 

2.  Surface- 


During  this  period,  a  weak  low  pressure  area 
originating  in  the  north  central  Pacific  deepened  and  moved 
to  the  northeast.  This  low  pressure  cell  eventually  reached 
the  British  Columbia  coast  at  the  end  of  the  period. 

A  closed  high  pressure  cell  was  located  in  the 
eastern  Pacific  off  the  northern  California  coast  ridging 
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Figure  6B 

1200  UTC  02  May  1990 
NOGAPS  Sea  Level  Pressure  Analysis 
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Figure  5D 

1200  UTC  03  May  1990 
NOGAPS  500  mb  Analysis 
heights  (solid,  ■> 
temperature  (dashed,  °C) 


Figure  6D 

1200  UTC  03  May  1990 
NOGAPS  Sea  Level  Pressure  Analysis 
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into  the  Pacific  Northwest.  The  high  pressure  pattern 
continually  weakened  while  gradually  moving  northeastward 
toward  the  Washington  coast. 

A  thermal  low  remained  over  Arizona  and  filled 
slightly  during  the  period.  This  low  had  an  associated 
inverted  trough  extending  over  the  California  coast.  (See 
Figures  6A  to  6D) . 

B.  MESOSCALE  FEATURES: 


Specific  mesoscale  features  of  interest  for  the  region 
during  the  period  include  the  land  and  sea  breeze,  the 
Catalina  eddy,  and  the  southerly  stratus  surge.  Satellite 
images  show  a  dominant  presence  of  stratus  along  the 
coastlines  during  this  time  period.  Because  of  the  close 
association  of  stratus  clouds  with  the  latter  two  features, 
this  study  concentrates  on  the  Catalina  eddy  and  the 
southerly  surge. 


1.  Catalina  Eddy- 

The  Catalina  eddy  is  a  feature  typically  occurring 
from  late  spring  to  early  fall  characterized  by  surface  wind 
cyclonic  circulation  in  the  vicinity  of  Santa  Catalina 
Island.  Figure  7,  from  Mass  et  al.  (1989),  is  a  Catalina 
eddy  composite  of  1200  UTC  surface  winds  (knots)  and  sea 
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Figure  7 

1200  UTC  Catalina  Eddy  and  climatology 
composite  of  surface  winds  (Knots)  and 
sea  level  pressure  (millibars). 

(Catalina  eddy  composite  based  on  50  events; 
climatology  composite  based  on  data  for  May 
through  September  1964-1982) 


(Mass  and  Albright,  1989) 
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Level  pressure  (mb)  compared  with  climatological  winds  and 


pressure . 

The  usual  coastal  southern  California  pattern  during 
this  period  is  characterized  by  westerly  to  northwesterly 
surface  winds  accompanied  by  morning  and  late  evening  fog 
and  stratus.  The  Catalina  eddy  will  often  form  when  coastal 
winds  shift  to  a  more  southerly  direction.  During  these 
eddy  conditions,  the  usual  dissipation  of  coastal  fog  and 
stratus  in  the  afternoon  will  often  not  happen. 

An  intensified  alongshore  pressure  gradient  created 
by  the  interaction  of  a  short  wave  trough  with  pre-existing 
troughs  is  credited  by  Mass  et  al.  (1989)  for  this  feature. 
Grandau  (1992)  found  this  eddy  to  extend  upward  to  about  920 
mb  for  this  particular  case.  These  eddy  events  can  vary 
from  a  few  hours  to  many  days  in  duration. 

2.  Southerly  Stratus  Surge- 

A  phenomena  along  the  west  coast  of  the  United 
States  that  can  often  be  associated  with  the  longer  duration 
Catalina  eddies  is  the  southerly  surge.  A  result  of  the 
alongshore  pressure  gradient  is  southerly  flow  developed 
within  about  100  km  (approximately  one  Rossby  radius)  of  the 
coastal  mountains.  As  the  southerly  winds  maintain  its 
flow,  the  cool  moist  marine  layer  is  deepened,  and  enhanced 
stratus  develops  near  the  coast.  The  deepened  marine  layer 
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will  often  result  in  improved  air  quality  ror  the  Los 
Ange  les  area . 

Figures  HA  through  8G  include  surface  wind 
observations  for  the  period  of  interest.  Gorki  1  l  (1991) 
described  a  weak  southerly  surge  for  his  study  covering  the 
same  time  period.  Figure  8F  shows  southerly  coastal  winds 
from  central  to  southern  California  at  2100  UTC  02  May  1990. 
The  2030  UTC  visible  satellite  image  (Figure  9)  reveals 
coastal  cloudiness  all  along  the  central  California  coast 
from  Monterey  to  San  Luis  Obispo  Bay.  Dorman  (1985) 
described  the  southerly  surge  as  a  result  of  coastal ly- 
trapped  gravity  currents,  while  Mass  et  al.  (1989)  stated 
that  it  results  from  a  coastal ly  trapped  two-layer  marine 
system. 
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Figure  9 

Visible  satellite  image 
(2030  UTC  02  May  90) 
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V.  CONTROL  MODEL  SIMULATIONS 


Stewart  (1992)  incorporated  the  Harshvardhan  et  al. 
(1987)  radiation  parameterization  into  the  NRL  regional 
weather  prediction  model  and  conducted  model  simulations 
integrated  for  36  hours  over  the  time  period  0000  UTC  02  May 
to  1200  UTC  03  May  1990.  Integration  over  the  same  area  and 
time  window  allows  for  consistent  comparisons  with  the 
control  model  of  Stewart  to  evaluate  the  impact  of  the  new 
stratus  parameterization. 

A.  CATALINA  EDDY- 

The  control  model  low-level  wind  fields  (Figures  8A  to 
8G)  indicated  clear  onset  of  a  Catalina  eddy  at  0900  UTC  02 
May  1990  with  south-southeasterly  winds  near  San  Diego  and 
clear  off-shore  flow  along  the  central  to  southern 
California  coast.  Further  off  the  coast,  general  northerly 
winds  were  observed  off  Monterey  Bay  to  west-northwesterly 
flow  toward  the  Mexico  border. 

Three  hours  later,  the  model  clearly  defined  the  eddy 
pattern  in  the  Los  Angeles  basin  with  the  vortex  centered 
between  the  Santa  Cruz  and  the  Santa  Catalina  Islands.  By 
1800  UTC,  the  pattern  had  become  disorganized  with  a 
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transition  or  southerly  winds  extending  along  the  coast  from 
San  Simeon  to  Santa  Barbara. 

At  2100  UTC,  the  winds  are  westerly  from  Pt .  Conception 
northward,  and  from  the  southwest  in  the  Los  Angeles  basin. 
By  0600  UTC,  another  eddy  onset  appeared  to  be  occurring, 
with  the  1200  UTC  wind  fields  again  giving  a  clear  eddy 
pattern . 

Comparison  of  model-predicted  wind  fields  and  a  limited 
number  of  land  and  ship  observations  showed  that  the 
Catalina  eddy  onset  appeared  to  be  reasonably  predicted. 
However,  the  model  seemed  to  show  the  dissociation  of  the 
eddy  sooner  than  the  observations  indicated.  The  1800  UTC 
reports  still  show  evidence  of  an  eddy  circulation,  while 
the  model  showed  mostly  onshore  winds  along  the  coast. 
Also,  observations  near  San  Diego  showed  the  winds 
maintaining  more  of  a  southerly  component  longer  into  the 
period  (through  0600  UTC  03  May) .  Because  of  the  lack  of 
data,  the  reformation  of  another  eddy  later  by  0600  UTC 
could  not  be  readily  verified. 
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B.  SOUTHERLY  SURGE - 


The  control  model  low-level  wind  fields  (Figures  8A  to 
8G)  showed  distinct  southerly  coastal  wind  flow  from  San 
Simeon  to  San  Diego  at  1800  UTC  02  May  1990.  By  2100  UTC . 
winds  from  San  Simeon  to  Santa  Barbara  shifted  to  a 
predominantly  western  direction.  From  0000  UTC  03  May  1990, 
winds  in  the  area  had  prominent  northerly  components. 

Observations  along  the  coast  showed  southerly  wind  flow 
from  Santa  Barbara  southward  at  1800  UTC  02  May  1990. 
Coastal  areas  to  the  north  however  did  not  have  southerly 
winds  until  2100  UTC  02  May  1990.  This  coastal  southerly 
wind  flow  continued  up  to  the  end  of  the  period. 

Coincident  with  the  southerly  surge  time  period, 
satellite  imagery  showed  persistent  cloudiness  along  the 
coastline  from  Monterey  Bay  to  San  Luis  Obispo  Bay.  The 
normal  late  day  coastal  fog  and  stratus  dissipation  occurred 
only  to  the  south.  Persistent  coastal  fog  and  stratus  with 
little  dissipation  is  expected  during  a  southerly  surge 
stage.  Figure  9  is  a  visible  satellite  image  highlighting 
the  coastal  cloudiness  for  the  region. 

Time-height  cross-sections  for  point  'A'  near  San  Luis 
Obispo  Bay  (35°N,  121  °W)  showed  model  southerly  winds  near 
the  surface  at  around  1800  UTC  02  May  time  period.  Figure 
10  shows  the  southerly  surface  flow  isolated  during  that 
time  window  (dashed  lines  indicate  northerly  winds  while 
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solid  lines  are  southerly) .  Also.  the  model  revealed  a 
local  maxima  in  humidity  in  the  lowest  200  meters  at  around 
1800  UTC  agreeing  with  the  expected  moist  marine  layer  that 
develops  from  the  southerly  surge.  Figure  11  highlights  the 
1  igh  moisture  near  the  surface  for  that  particular  time 
per lod . 
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VI.  TEST  MODEL  RUN  EVALUATIONS 

The  test  model.  which  utilizes  the  Sundqvist  et  al. 
explicit  non-convect 1 ve  cloud  water  scheme,  was  compared 
with  the  control  model,  Stewart's  (1992)  incorporation  of 
the  Slingo  and  Ritter  (1985)  scheme  into  the  NRL  mesoscale 
model.  Specific  features  compared  were  mesoscale  wind 
strucuture.  cloud  fraction,  longwave  radiation.  and 
temperature  for  the  period  0000  UTC  02  May  1990  to  1200  UTC 
03  May  1990.  In  addition,  the  depiction  of  cloud  structure 
by  cloud  liquid  water  is  evaluated  and  compared  with  that  of 
cloud  fraction.  The  evaluations  focused  on  the  cloud, 
temperature,  radiation,  and  moisture  structure  within  the 
boundary  layer. 

A.  MESOSCALE  WIND  FEATURES- 

Low-level  wind  fields  defining  the  control  model 
simulations  of  the  Catalina  eddy  and  southerly  surge  are 
described  in  Section  V.  Comparisons  of  test  model  to 
control  model  low-level  winds  revealed  differences  of  less 
than  0.5  m/s  between  the  fields.  This  may  be  explained  by 
the  incorporation  of  the  cloud  liquid  water  scheme 
predominantly  affecting  smaller  scale  features.  The  scheme 
also  affects  thermodynamic  processes  more  than  dynamic  ones. 
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The  j  'at  all  na  <-ddy  arid  -southerly  surge  occur  near  the 
coastlines  where  the  cloud  cover  test  model  simulations 
showed  little  change  when  compared  to  the  control  model. 

Because  of  the  similarities  between  these  mesoscaie  wind 
fields,  the  following  comparisons  will  *ocus  on  cloud 
fraction,  temperature,  and  longwave  radiation  outputs. 

B.  CLOUD  FRACTION- 

Comparison  of  the  control  model  to  test  model  cloud 
fractions  showed  the  test  model  to  be  more  deficient  than 
the  control  case  in  forecasting  cloud  cover  over  the  region 
as  verified  against  satellite  observations.  In  fact  the 
test  model  for  2100  UTC  02  May  simulated  no  clouds  at  the 
1000  mb  leve 1 . 

Since  boundary  layer  depth  has  a  direct  effect  on  the 
relative  humidity  threshold  values  (RHC)  in  the  test  model, 
the  location  of  the  boundary  layer  top  appears  to  be  the 
major  reason  for  the  test  model's  weakness  in  depicting 
cloud  regions.  Figure  12  illustrates  the  effects  of  a 
shallow  versus  a  deep  boundary  layer  on  cloud  formation. 
Because  of  the  assumption  of  a  linear  decrease  of  RHC  with 
height  to  the  top  of  the  boundary  layer,  a  shallow  boundary 
layer  results  in  lower  RHc  values  within  the  layer  making  it 
easier  to  produce  clouds  at  those  lower  levels.  On  the 
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Figure  12 

Effect  of  Boundary  Layer  Depth 
on  cloud  formation. 
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other  hand,  a  deep  layer  results  in  more  difficulty 
producing  clouds. 

Also  shown  in  Figure  12  is  the  height  at  which  Slingo 
and  Ritter's  lower  RHc  profile  coincides  with  the  Sundqvist 
et  al.  profile  (over  water) .  This  occurs  when  the  top  of 
the  boundary  layer  is  at  a  sigma  level  of  about  0.96. 
Therefore  a  boundary  layer  depth  less  than  approximately  400 
meters  is  required  for  the  test  model  to  produce  clouds 
before  the  control  model.  Examination  of  boundary  layer 
structure  over  the  ocean  for  this  case  generally  showed 
depths  greater  than  400  meters.  Hence,  one  would  not  expect 
increased  low-level  cloud  formation  in  the  test  model. 

Stewart  (1992)  described  the  inability  to  diagnose 
certain  cloudy  regions  as  one  of  the  weaknesses  of  the 
control  model,  resulting  in  his  modification  of  the  critical 
relative  humidity  profile  at  the  low  levels  (the  dashed  line 
of  Figure  4) .  Stewart's  profile  allows  more  cloud  formation 
at  the  lower  levels.  Thus,  henceforth  the  control  model 
will  be  compared  with  two  additional  model  runs  —  one  using 
the  Sundqvist  et  al.  explicit  cloud  liquid  water  scheme  but 
retaining  the  Slingo  critical  relative  humidity  profile 
(test  model  ’A’),  and  one  using  the  Sundqvist  et  al.  scheme 
incorporating  Stewart's  proposed  modified  profile  (test 
mode  1  ' B ' )  . 
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i:  Cloud  Fraction  (Control  Model  vs  Test  Model _ '  A  1 j - 


Comparison  of  1000  mb  cloud  fractions  for  2100  UTC 
02  May  for  this  case  (Figure  13)  shows  that  test  model  'A' 
significantly  extended  the  cloud  region  westward  and 
southward  from  approximately  32°N.  125°  W.  However.  there 
was  minimal  effect  toward  the  California  coastline.  The 
differences  in  the  cloud  fractions  are  due  to  differences  in 
relative  humidity.  The  inclusion  of  the  Sundqvist  et  al. 
(1989)  cloud  liquid  water  scheme  causes  a  temperature 
decrease  in  regions  of  cloud  liquid  water.  The  temperature 
drop  in  these  regions  consequently  results  in  higher  RH 
values.  Temperature  effects  are  discussed  further  in  Part  D 
of  this  Section. 

Figure  14  compares  RH  values  at  1000  mb.  The 
increased  relative  humidities  toward  the  southwest  quadrant 
account  for  the  increased  cloud  cover  in  that  region.  The 
coastal  relative  humidities  on  the  other  hand  remain  fairly 
constant  resulting  in  minimal  cloud  cover  increase  near  the 
coast . 
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Figure  13 

1000  ab  cloud  fraction 
Control  model  -  top 
Test  model  ’A’  -  bottom 


2100  UTC  02  May  90 


iZO  113 

(Contour  interval  = 5 %) 


Figure  14 

1000  mb  relative  humidity, 
(control  model  —  solid) 
(test  model  *A*  —  dashed) 
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2.  Cloud  Fractions  (Control  Model  vs  Test  Model  ' B ' ) — 


Test  model  'B'.  in  addition  to  extending  the  low- 


level  cloud 

region  westward 

s imi 1 ar 

to  case  '  A  '  , 

a  Iso 

realistical ly 

extended 

the 

c  loud 

region 

toward 

the 

California  coast  (Figure 

15)  . 

Since 

Stewart 

mod if  led 

RHC 

values  to  help  produce  more  model-simulated  clouds  by  using 
satellite  imagery  and  observations,  this  extension  of  clouds 
along  the  coastline  can  be  expected. 

Satellite  imagery  (2030  UTC  visible  image.  Figure  9) 
showed  a  low-level  tongue  of  coastal  cloudiness  from 
Monterey  Bay  to  San  Luis  Obispo  Bay  that  the  model  did  not 
simulate.  In  addition,  the  model  predicts  low-level  cloud 
cover  out  about  100  run  seaward  from  the  coast  (34°  N,  123°  W) 
where  generally  clear  skies  were  observed.  The  sea  surface 
temperature  (SST)  fields  used  in  these  simulations  were  from 
0000  UTC  02  May  1990  and  kept  constant  throughout  the  time 
period.  The  constant  SST  fields  may  be  part  of  the  reason 
for  the  model  cloud  fraction  weaknesses.  SST  values  warmer 
than  or  approximately  equal  to  the  surface  air  temperature 
may  result  in  more  low-level  clouds  due  to  increased 
boundary  layer  flux  convergence  while  significantly  cooler 
model  SST  values  could  inhibit  cloud  formation. 

As  discussed  in  Section  V,  the  control  model's  (and 
subsequently  the  test  model's)  simulation  of  low-level  winds 
resulted  in  a  much  shorter  duration  southerly  surge  event  as 
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compared  with  wind  observations  along  the  California  coast. 
Since  the  southerly  surge  often  produces  persistent  coastal 
cloudiness,  the  weaker  model  simulation  of  southerly  winds 
may  be  another  cause  of  the  model's  weakness  in  producing 
clouds  up  to  the  coastline. 

The  increased  mode  1 -s imu 1 ated  cloud  cover  should 
affect  the  radiation  and  temperature  fields  at  and  below  the 
cloud  levels.  This  is  evaluated  in  the  following  sections. 

C.  LONGWAVE  RAD  I  ATI 0N- 

Figure  16  depicts  an  idealized  longwave  radiation 
profile  in  relation  to  a  stratocumu lus  cloud  (Stull  1988). 
A  maximum  longwave  cooling  rate  is  expected  at  the  top  of 
the  cloud  boundary  with  minimal  cooling  above  the  cloud. 
Also  weaker  longwave  heating  may  be  expected  at  the  cloud 
base . 

Longwave  radiation  along  with  cloud  fraction  and  cloud 
liquid  water  profiles  are  shown  in  Figure  17  for  Point  'B' 
(33“  N,  122“  W)  for  test  model  'A'  at  2100  UTC  02  May.  Cloud 
tops  occur  at  about  990  mb  for  both  cloud  fraction  and  cloud 
liquid  water  representations.  Also,  cloud  fraction  and 
cloud  liquid  water  profiles  depict  the  cloud  base  at  about 
1005  mb.  Maximum  longwave  cooling  of  approximately  78°C/day 
occurs  within  the  cloud  at  about  1000  mb  with  a  heating  rate 
of  about  17°C/day  at  the  cloud  bottom.  Evaluating  this 
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Figure  10 

Idealized  longwave  radiation  cooling  rate 
profile  In  a  stratocumulus  deck. 

(Stull  1988) 
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Figure  17 

Longwave  radiation  cooling  profile 
at  Point  ’B* 

(Test  model  ’A*,  2100  UTC  02  May  90) 
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model  radiation  profile  with  the  ideal  case  represented 
Figure  16  shows  that  the  model's  maximum  longwave  radiat 
cooling  is  t epresented  favorably  though  the  model  ioca 
the  cooling  maximum  at  a  lower  level  within  the  clo' 
Incorporation  of  the  cloud  liquid  water  scheme  provide 
direct  interaction  of  the  cloud  with  the  environmental 
temperature  and  moisture  profiles.  This  interaction  is 
evident  in  the  equations  (3.1.  3.6.  and  3.7)  as  given  by 
Sundqvist.  Thus  changes  in  cloud  liquid  water  as  computed 
at  every  griapoint  and  for  each  iteration  in  the  model 
directly  impact  the  temperature  and  moisture  of  the 
atmosphere . 

However,  these  changes  in  cloud  liquid  water  are  not  as 
readily  apparent  when  comparing  cloud  liquid  water  to  cloud 
fraction.  This  is  because  cloud  fraction  is  computed  along 
with  radiation  parameters  every  half  hour  (60  iterations) . 
Thus  the  indirect  effect  of  cloud  liquid  water  on  cloud 
fraction  through  changes  in  temperature  and  moisture 
profiles  is  not  readily  visible  in  comparisons.  Hence, 
cloud  liquid  water  provides  a  more  physically  realistic 
determination  of  cloud  structure  and  subsequent  radiation 
profiles,  though  for  comparisons  of  model  simulations  there 
is  a  more  direct.  one-to-one  relationship  betwen  cloud 
fraction  and  radiation  than  cloud  liquid  water  and 
radiation . 
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D.  TEMPERATURE - 


Temperature  field  comparisons  revealed  significant 
differences  between  the  control  model  and  test  models  A' 
and  1 B ' .  The  temperature  differences  occurred  within  the 
mode 1-s imu 1 ated  cloud  regions  for  both  test  models. 

1.  Temperature  (Control  Model  vs  Test  Model  'A')- 

Figure  18  shows  low-level  cross-section  temperature 
comparisons  of  the  control  model  and  test  model  'A'  for  a 
predominantly  cloud-covered  region  from  Point  ’ C'  to  ‘D’  at 
2100  UTC  02  May.  Temperatures  were  as  much  as  2°C  lower 
for  test  model  'A'  versus  the  control  model  in  areas  within 
and  below  the  dense  clouds.  Figure  19  for  Point  'B'  shows 
that  temperature  deviations  occurred  from  the  surface  up  to 
about  960  mb  with  the  maximum  difference  occurring  within 
the  cloud  boundaries. 

2.  Temperature  (Control  Model  vs  Test  Model  ‘ B ' ) — 

Low-level  cross-section  temperature  comparison  for 
the  same  region  described  above  is  shown  in  Figure  20. 
Similar  to  test  model  'A'.  temperatures  in  test  model  'B' 
were  lower  by  as  much  as  4° C  within  cloud  regions  as 
compared  to  the  control  model.  The  temperature  field 
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(Contour  interva 1 = 1 °K ) 

Figure  18 

Temperature  cross-section  from  Point  *C’  to  ’D' 
(2100  UTC  02  May  90) 

Control  model  -  solid 
Test  model  ’A*  -  dashed 


Figure  19 

Temperature  profile  at  Point  'B' 
(2100  UTC  02  May  90) 
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(Contour  interval=l°K 


Figure  20 

Temperature  cross-section  from  Point  ’ 
(2100  UTC  02  May  90) 
Control  model  -  solid 
Test  model  ’  B’  -  dashed 
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to  'D 


difference  for  this  case  appeared  to  be  larger  due  to 


thicker  cloud  depictions  by  the  model. 

Also  comparable  to  test  model  'A',  Figure  21  for 
Point  ' B'  shows  temperature  profile  deviation  occurring 
within  the  same  cloud  boundary  levels. 

E.  CLOUD  LIQUID  WATER- 

Generally,  low-level  cloud  liquid  water  regions  closely 
follow  the  cloud  fraction  regions  as  seen  in  Figure  22. 
Note  the  rather  smooth  nature  of  cloud  fraction  (as 
dependent  on  RH)  in  contrast  to  the  somewhat  noisy  structure 
of  cloud  liquid  water. 

An  area  of  interest  due  to  the  lack  of  clouds  at  1000  mb 
is  located  in  the  region's  southwest  quadrant  (Figure  23). 
A  comparison  between  the  test  model  run  'B'  1000  mb  cloud 
fraction  and  cloud  liquid  water  fields  for  2100  UTC  02  May 
show  little  or  no  cloud  liquid  water  in  that  quadrant  at 
1000  mb.  A  cross  section  from  point  'C'  to  'D'  (Figure  24) 
shows  cloud  liquid  water  generally  throughout  that  region 
with  relatively  dry  regions  at  the  1000  mb  level.  This 
layered  cloud  structure  is  often  observed  in  the  marine 
boundary  layer.  Use  of  the  explicit  cloud  liquid  water 
scheme  allows  for  a  realistic  depiction  of  this  layered 
structure  that  a  simple  cloud  fraction  scheme  could  not 
simulate.  Note  how  the  cloud  liquid  water  cross  section 
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Figure  21 

Temperature  profile  at  Point  '  B* 
(2100  UTC  02  May  90) 
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Cloud  Liquid  Water 
contour  interval=0 


Figure  23 

1000  mb  cloud  fraction  (top) 

1000  mb  cloud  liquid  water  (bottom) 
(Test  model  *  B'f  2100  UTC  02  May  90) 
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contour  intervals: 

Cloud  liquid  water=0.1  g/kg 
Cloud  fractions. 1 


Horizontal  distance  *  870.8  km 


Figure  24 

Cloud  liquid  water  (solid)  and 
cloud  fraction  (dashed)  cross-sections 
from  Point  'C'  to  ’ D’  . 

(Test  model  ’B’,  2100  UTC  02  May  90) 


68 


also  provides  a  more  realistic  horizontal  ana  vertical  cloud 
texture  appearance  than  the  cloud  fraction  contours  as 
compared  in  Figure  24.  The  cloud  fraction  is  simply  a 
representation  of  RH  while  cloud  liquid  water  includes  more 
of  the  physics  needed  in  defining  clouds  and  the  boundary 
1 ayer . 


69 


VII.  '-'QNCLUS I ONS  AND  RECOMMENDATIONS 


A.  CONCLUS I ONS- 

This  study  has  demonstrated  that  the  incorporation  of 
the  Sundqvist  et  al.  (1989)  explicit  non-convect 1 ve  cloud 
liquid  water  scheme  into  the  NRL  mesoscale  model  provides 
improvements  in  simulating  longwave  radiation,  temperature, 
and  cloud  structure  features  over  the  U.S.  West  Coast 
Pacific  region. 

One  limitation  that  was  found  in  the  incorporation  of 
the  Sundqvist  scheme  for  this  particular  case  was  its 
extreme  weakness  in  accurately  diagnosing  low-level  clouds. 
The  main  reason  for  this  weakness  was  Sundqvist 's 
representation  of  relative  humidity  threshold  values  for 
cloud  formation.  Test  model  evaluations  were  therefore 
conducted  using  Sundqvist 's  cloud  liquid  water  scheme  but 
retaining  Slingo's  critical  relative  humidity  profile  for 
cloud  formation  (test  model  'A').  In  addition,  test  model 
' B ' ,  which  used  a  modified  critical  relative  humidity 
profile  by  Stewart  (1992)  was  also  evaluated.  This  modified 
profile  was  proposed  to  help  compensate  for  the  control 
model's  weakness  in  diagnosing  low-level  clouds  over  certain 
areas  in  the  region.  Test  models  'A'  and  'B'  were  able  to 
provide  a  closer  cloud  cover  picture  of  the  region  as 
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compared  to  satellite  imagery.  with  test  model  B ‘  showing 
more  realistic  cloud  cover  near  the  coast. 

Longwave  radiation  profiles  in  cloud  regions  showed 
realistic  and  consistent  cooling  and  warming  rates  as 
compared  to  idealized  radiation  profiles.  Temperature 
effects  due  to  model-simulated  cloud  cover  also  compared 
f  avorab  ly . 

A  significant  improvement  in  the  NRL  model  was  the 
realistic  horizontal  and  vertical  cloud  structure  that  was 
represented  by  the  cloud  liquid  water  fields.  Control  model 
cloud  depiction  was  only  through  cloud  fraction,  a 
representation  of  relative  humidity.  The  introduction  of 
cloud  liquid  water  as  a  prognostic  variable  takes  into 
account  more  of  the  physics  involved  in  better  defining 
cloud  structure  and  the  marine  boundary  layer.  Cloud  liquid 
water  used  as  a  portrayal  of  cloud  cover  provides  a  more 
physically  realistic  texture  and  layered  structure 
associated  with  non-convect 1 ve  clouds. 
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B.  RECOMMENDATI ONS- 


Stewart ' s  (1992)  study  provided  a  modified  RHC  profile 
to  help  provide  a  closer  cloud  depiction  as  compared  to 
satellite  imagery.  Using  this  modified  profile,  better 
cloud  cover  representation  was  achieved  but  continued  fine- 
tuning  of  the  RHC  values  for  a  variety  of  cases  of  differing 
synoptic  flow  may  be  necessary  for  a  more  accurate  model 
cloud  cover  picture.  An  experiment  utilizing  a  good 
dispersal  of  surrace.  ship,  and  buoy  station  observations 
with  soundings  along  with  satellite  observations  could 
provide  a  closer  real  world  RHC  profile  depiction. 

Section  VI  described  a  situation  where  the  1000  mb  cloud 
liquid  water  depiction  showed  very  little  or  no  clouds  while 
satellite  imagery  clearly  indicated  low-level  cloud  cover 
over  the  area.  A  cross-section  view  of  the  test  model  cloud 
liquid  water  fields  revealed  a  layered  cloud  structure  with 
a  dry  area  at  the  1000  mb  level.  Development  of  a  three- 
dimensional  display  capability  would  help  provide  an  easier 
way  of  visualizing  a  complex  cloud  structure  from  cloud 
liquid  water.  Computation  of  cloud  fraction  at  every 
iteration  to  correspond  to  cloud  liquid  water  fields  would 
aid  in  visual  interpretation  of  clouds.  In  addition.  other 
fields  in  the  model  may  be  easier  visualized  three- 
dimensiona 1 ly . 
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This  study  incorporated  the  Sunaqvist  et  al.  >1989)  non- 


convect  1  ve  cloud  scheme  into  the  NRL  model  but  did  not 
include  Sundgvist's  scheme  tor  convective  clouds.  Further 
study  incorporating  the  convective  cloud  scheme  is 
recommended  as  a  follow  on. 

Finally,  an  evaluation  of  the  realistic  cloud  structure 
by  cloud  liquid  water  in  regions  of  available  data  could 
verify  the  accuracy  of  the  cloud  portrayals  and  provide  a 
better  degree  of  confidence  in  the  overall  model  outputs. 
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